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ABSTRACT 

Context. X-ray flashes (XRFs) are a class of gamma-ray bursts (GRBs) with the peak energy of the time-integrated v spectrum, 
£p, typically below 30 keV, whereas classical GRBs have Ep of a few hundreds keV. Apart from E^, and the systematically lower 
luminosity, the properties of XRFs, such as the duration or the spectral indices, are typical of the classical GRBs. Yet, the nature of 
XRFs and the differences from that of GRBs are not understood. In addition, there is no consensus on the interpretation of the shallow 
decay phase observed in most X-ray afterglows of both XRFs and GRBs. 

Aims. We examine in detail the case of XRF 080330 discovered by Swift at the redshift of 1.51. This burst is representative of the 
XRF class and exhibits an X-ray shallow decay. The rich and broadband (from NIR to UV) photometric data set we collected across 
this phase makes it an ideal candidate to test the off-axis jet interpretation proposed to explain both the softness of XRFs and the 
shallow decay phase. 

Methods. We present prompt y-ray, early and late NIR/visible/UV and X-ray observations of the XRF 080330. We derive a spectral 
energy distribution from NIR to X-ray bands across the shallow/plateau phase and we describe the temporal evolution of the multi- 
wavelength afterglow within the context of the standard afterglow model. 

Results. The multi-wavelength evolution of the afterglow is achromatic from ~10^ s out to ~8x 10"* s. The energy spectrum from NIR 
to X-ray is nicely fitted with a simple power-law, F,, oc v^^"", with /?ox = 0.79 + 0.01 and negligible rest-frame dust extinction. The 
light curve can be modelled either by a piecewise power-law or by the combination of a smoothly broken power law with an initial 
rise up to ~600 s, a plateau lasting up to ~2 ks, followed by a gradual steepening to a power-law decay index of ~2 out to 82 ks. At 
this point, there appears a bump modelled with a second component, while the corresponding optical energy spectrum, Fy oc v^^" , 
reddens by A/3o = 0.26 ± 0.06. 

Conclusions. A single-component jet viewed off-axis explains the light curve of XRF 080330, the late time reddening as due to the 
reverse shock of an energy injection episode, and its being an XRF. Other possibilities, such as the optical rise marking the pre- 
deceleration of the fireball within a wind environment, cannot be definitely excluded, but seem somewhat contrived. We rule out the 
dust decreasing column density swept up by the fireball as the explanation of the rise of the afterglow. 

Key words, gamma rays: bursts; X-rays: individual (XRF 080330) 



1. Introduction 

Time-integrated photon spectra of long gamma-ray bursts 
(GRBs) can be adequately fitted with a smoothly broken power 
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law jBand et al.lll993h . whose low-energy and high-energy in- 
dices, aB and 13b, have median values of -1 and -2.3, respec- 
tively dPreece et al.ll200(A iKaneko et alj|200 6l). The correspond- 
ing vFy spectrum peaks at Ep, the so-called peak energy, whose 
rest frame value is found to correlate with other relevant ob- 
served intrinsic properties, such as the isotropic -equivalent ra- 
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dialed y-ray energy, Ejso ( Amati et alj 20021) , or its coUimation- 
corrected value, Ey ('Ghirla nda eTall 12004'). In the BATSE cat- 
alogue, the Ep distributi on clusters around 300 keV with a 
-100 keV width ( Kaneko etai]|2006l) . 

When observations of GRBs in softer energy bands than 
BATSE became available thanks to BeppoSAX and HETE- 
2, a new class of soft GRBs with Ej,< 3 keV, so 
name d X-ray flashes (XRF s), was soon discovered dHeise et alj 
120011: iBarraud et all l2003l) . GRBs with intermediate softness, 
called X-ray rich (X RR) bursts, were als o observed, with 
30 keV ^ ^100 keV dSakamotoet al.l2005h . These soft GRBs 
share the same temporal and spectral properties, aside from the 
systemat ically lower E^, with the classical GRBs both for the 
prompt dFroriteraetal J 120001: iBarraud et all 120031: lAmati et al " 



2005 



2004 ) and, partially, the afterglow emission ("Sakamo to et al 



ID' Alessio et al.l 2006: Mangano et al. 2007). M oreover 

they were found to obey the Ep-Ei^o correlation (Amati et alj 
|2007|) discovered for classical GRBs, extending it all the way 
down to iSp values o f a few keV and forming a continuum 
dSakamoto et alj|2005l) . Like classical GRBs, also XRFs have 
been found to b e associated with SNe (Campana et al. 2006: 
iPian et al] l2006h and therefore connected with the collapse of 
massive stars. The comparison also holds for the cases in which 
apparently no associated SN was found both for classical GRBs 
(e.g. DeUa Valle et al. 2006, Fynbo et al. 2006, Gal-Yam et al. 
2006) and for XRFs (Levan et al. 2005). 

A number of different models have been proposed in the 
literature to explain the nature of XRRs and XRFs (e.g., see 
the review of Zhang 2007): i) standard GRBs viewed well off 
the axis of the jet, thus explaining the softn ess as due to a 
larger viewing angle and a lowe r Doppler factor dYamazaki et alj 
120021: iGranot et alj|2002l |2005|) : ii) two coaxial jets with differ- 
ent opening angles (wide a nd narrow), 0^, > 0^, and viewed at 
an angle 0^, 0,, < 0y < 0„ ( Peng et al.ll2005h : iii) the "dirty fire- 
ball" model characterised by a small value of the bulk Lorentz 
factor due to a rela tively high baryon loading of the fireball 
dPermer et alj l2000l) : iv) distribution of hi gh Lorentz factors 
with low contrast of the colliding shells (" Mochkovitch et al] 
12004 . In the off-axis interpretation, a number of different mod- 
els of the structure and opening angle of the jet have been pro- 
posed (e.g. Granot et al. 2 005; Donaghy 2006 ). 

The advent of Swift dGehrels et alj|2004l) has made it pos- 
sible to collect a large sample of early X-ray afterglow light 
curves of GRBs. Concerning the XRRs and XRFs, the 15- 
150 keV energy band of the Swift Burst Alert Telescope (BAT; 
Barthelmy et al. 2005) and its relatively large effective area 
still allow to detect them, although those with Ep of a few keV 
are disfavoured with respect to BeppoSAX and HETE-2 instru- 
ments (Sakamoto et al. 2008). Thanks to Swift it is possible to 
study the early X-ray afterglow properties of these soft events. 
Like hard GRBs, als o XRFs occasionally ex hibit X-ray flares 
dRomano et al.ll2006h . ISakamoto et all d2008h analysed a sam- 
ple of XRFs, XRRs and classical GRBs detected with Swift and 
found some evidence for an average X-ray afterglow luminos- 
ity of XRFs being roughly half that of classical GRBs and some 
differences between the average X-ray afterglow light curves. 

An unexpected discovery of Swift is the shallow decay phase 
experienced by most of X-ray afterglows bet ween a few 10^ 
up to 10^-10'* s af t er the trigger time (.Tagliaferri et al.l l2005t 
iNousek et al] 120061; Izhang et al.l l2006l) . Several interpretations 
have been put forward (e.g. see Ghisellini et al. 2009 for a 
brief review). Among them, some invoke continuous energy 
injection to the fireball shock front t hrough refreshed shocks 
dNousek et alJl2006l : IZhang et ani2006i) . depending on the pro- 



genitor period of activity: a long- or short-lived powering mech- 
anism, either in the form of a prolonged, continuous energy re- 
lease (L(f) ocr*), or via discrete shells whose F distribution is 
a steep power-law. For instance, in the cases of GRB 050801 
and GRB 0701 10 a newly born millisecond magnetar was sug- 
gested to power the flat decay observed in the op tical and X-ray 
bands (Pe Pasqual e et al.ll2007HTroja et alJl2007i) . Alternatively, 
geometrical models interpret the shallow decay as the delayed 
onset of the a fterglow observed from viewing ang l es outside the 
edge of a iet dGranot et al.ll2002HSalmonsonl2003trGranoll2005l: 



lEichler & Grano j 120061) . Other models invoke two-component 
jets viewed off the axis of the naiTow component, also inv oked to 
expla in the late time observations of GRB 030329 (Berger etjlj 
|2003|) . In particular, this model would explain the initial flat de- 
cay observed in XRF 030723, dominated by the wide compo- 
nent, followed by a late rebrightening peaking at ~16 days and 
interpreted as due to th e deceleration and late ral expansion o f 
the narrow component dHuang et alj[2004; .Butler et a n i2005h . 
although alternative explanations for this in terms of a SN h ave 
also been proposed dFvnbo et alJl2004l ; lTominaga et al.ll2004l) . 

Other models explain the shallow decay as due to a temporal 
evolu t ion of the fireball micro-physical parameters dioka et al.l 
120061 : iGranot et alJ '2006'): scattering by dust located in the cir- 
cumburst medium (Shao & Dai 2007); "late prompt" activity of 
the inner engine, keeping up a prolonged emission of progres- 
sively lower power and Lorentz factor shel ls, which radiate at the 
same distance as for the prompt emission dGhisellini et al.ll2b07l 
,20091) ; a dominating reverse shock in the X-ray band propagating 
through late shells with sma l l Lorentz factors ( Genet et al. 20071; 
lUhm & Beloborodovll2007h . lYamazakil d2009l) suggests that the 
plateau and the following standard decay phases are an artifact of 
the choice of fo, provided that the engine activity begins before 
the trigger time by ~10"'-10'* s. 

XRF 080330 was promptly discovered by the Swift-BAT 
and automatically pointed with the X-Ray Telescope (XRT; 
Burrows et al. 2005) and the Ultraviolet/Optical Telescope 
(UVOT; Roming et al. 2005) as shown in Figure |3] In this work 
we present a detailed analysis of the Swift data, from the prompt 
y-ray emission to the X-ray and optical afterglow and combine 
it with the large multi-filter data set collected from the ground, 
encompassing a broad band, from NIR to UV wavelengths, and 
spanning from one minute out to ~3 days post burst. The main 
properties exhibited by XRF 080330 are the rise of the optical af- 
terglow up to -300 s, followed by a shallow decay also present in 
the X-ray, after which it gradually steepens, and either a possible 
late time (~ 10^ s) brightening (Fig.|5]l or a sharp break (Fig.|4]i. 
The richness of the multi-wavelength data collected throughout 
the rise-flat top-steep decay allows us to constrain the broadband 
energy spectrum of the shallow decay phase as well as its spec- 
tral evolution. Moreover, it is possible to constrain the optical 
flux extinction due to dust along the line of sight and, in par- 
ticular, near the progenitor. This GRB is a good benchmark for 
the proposed models of XRFs sources and of their link with the 
classical GRBs through the common properties, such as the flat 
decay phase. 

The paper is organised as follows: Sects. |2] and |3] report 
the observations, data reduction and analysis, respectively. We 
report our multi-wavelength combined analysis in Sect. |4] In 
Sect.|5]we discuss our results in the light of the models proposed 
in the literature and Sect. |6]reports our conclusions. 

Throughout the paper, times are given relative to the BAT 
trigger time. The convention F(v, t) oc v"^ r" is followed, where 
the spectral index /? is related to the photon index F = [i + \. 
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We adopted the standard cosmology: //q = 70kms ' Mpc \ 
Qa = 0.7, Qm = 0.3. 

All the quoted errors are given at 90% confidence level for 
one interesting parameter {Ay^ - 2.706), unless stated other- 
wise. 



tinction in each filter was estimated through the NASA/IPAC 
Extragalactic Database extinction calculator: Auvw\ - 0.120, 
Au = 0.090, Ab = 0.071, = 0.064, Ay = 0.055, A, = 0.047, 
Ar = 0.044, A/ = 0.032, A, = 0.035, A, = 0.022, Ay = 0.015, 
Ah = 0.010, Ax =0.006. 



2. Observations 

XRF 080330 triggered the Swift-BAT on 2008 March 30 at 
03:41:16 UT. The y-ray prompt emission in the 15-150 keV en- 
ergy band cons isted of a m ultiple-peak structure with a duration 
of about 60 s jMao et alj |2008a). An uncatalogued, bright and 
fading X-ray source was promptly identified by XRT. From the 
initial 100-s finding chart taken with the UVOT telescope in the 
White filter from 82 s the optical counterpart was initially lo- 
calised at RA = ll'' 17™ 04!51, Dec. = -h30° 37' 22'.'1 (J2000), 
with an error radius of I'.'O (Icr; Mao et al. 2008a). During the 
observations, the Swift star trackers failed to maintain a proper 
lock resulting in a drift which affected the observations and ac- 
curacy of early reports. We finally refined the position from the 
UVOT field match to the USNO-Bl catalogue: RA = ll'' 17™ 
04=:52, Dec. = -h30° 37' 23'.'5 (J2000), with an error radius of 
0'.'3 (Icr; Mao et al. 2008b), consistent with the position derived 
from ground telescopes (e.g., PAIRITEL, Bloom & StaiT 2008). 

The Telescopes a Action Rapide pour les Objets Transitoires 
(TAROT; Klotz et al. 2008c) began observing at 20.4 s (4.5 s 
after the notice) and discovered independen tly the optical coun - 
terpart during the rise with /?~16.8 at 300 s dKlotz e t al."2008a'). 
TAROT went on observing until the dawn at 1 .4 ks ( Klotz et al, 
I2008bl) . 

The Rapid Eye MounQ (REM; Zerbi et al. 2001) telescope 
reacted promptly and began ob serving at 55 s and detected the 
optical afterglow in R band dD'Avanzo et al.l 12008). The opti- 
cal counterpart was promptly d etected also by other r obotic tele- 
scope s, such as R OTSE-IIIb dSchaefer et al.l 120081; lYuan et alj 
12008 1). PROMPT ( ISchubeletalJ l2008f) and RAPTOR; the lat- 
ter in particular observed a ~10-s long optical flash of R - 
17.46 ± 0.22 at 60 s contemporaneous with the last y-ray pulse 
(IWren et alj|2008) . 

The Liverpool Telescope (LT) began observing at 181 s. 
The optical afterglow w as au tomatically identified by the 
LT-TRAP GRB pipeline dGuid orzi et al. 2006) with r'~17.3 
dGomboc et alj|2008al) . thus triggering the multi-colour imaging 
observing mode in the g'r'i' filters which lasted up to the dawn 
at 4.9 ks. 

The FauLkes Telescope North (FTN) observations of 
XRF 080330 were caiTied out from 8.4 to 9.1 hr and again from 
31.8 to 33.9 hr with deep r' and /' filter exposures. 

The Gamma-Ray Burst Optical and Near-Infrared Detector 
(GROND; Greiner et al. 2008) started simultaneous observations 
in g'r'i'z'JHK filters of the field of GRB 080330 at 3.1 min- 
utes and detected the afterglow with J - 15.92 + 0.04 and 
H = 15.46 ±0.11 from the first 240-s of effective exposure 
dClemens et al ] l2008l) . 

A spectrum of XRF 080330 was acquired at 46 minutes with 
the Nordic Optical Telescope (NOT). The identification of ab- 
sorption features all owed to measure the redshift, which turned 
out to be z = 1.51 dMalesani et al.ll2008l) . This was soon con- 
firmed by the spectra ta ken with the Hobby-Eberly Telescope 
dCucchiara & Foxll2008l) . 

The Galactic redd ening along the line of sight to the GRB 
is Eb-v - 0.017 (.Schlegel et aljfl998i) . The corresponding ex- 

' |http : / /www ■ r em ■ inaf . i t/l 



3. Data reduction and analysis 

3.1. Gamma-ray data 

The BAT data were processed with the heasoft package (v.6.4) 
adopting the ground-re fined coordinates provided by the BAT 
team (M arkwardt et al.li2008.) . The BAT detector quality map 
was obtained by processing the nearest-in-time enable/disable 
map of the detectors. 

The top panel of Figure [T] shows the 15-150 keV mask- 
weighted light curve of XRF 080330 as recorded by BAT, ex- 
pressed as counts per second per fully illuminated detector for 
an equivalent on-axis source. The solid line displayed in Fig. [1] 
coiTesponds to the result of fitting the profile from -1.2 to 100 s 
with a combination of four pulses (Markwardt et al. 2008) as 
modelled by Norris et al. (2005; hereafter N05 model). Table [1] 
reports the corresponding derived parameters: fp (peak time), A 
(15-150 keV peak flux), (rise time), (decay time), w (pulse 
width), k (pulse asymmetry) and the model fluence in the 15- 
150 keV band. The goodness of the fit is ^-^dof = 375/379. 
Parameter uncertainties were derived by propagation starting 
from the best-fit parameters and taking into account their covari- 
ance. We tried to apply the same analysis to the light curves of 
the resolved energy channels to investigate temporal lags and, 
more generally, the dependence of the parameters on energy; 
however, because of the faintness and softness of the signal, we 
could not constrain the parameters in a useful way. 

In addition, the energy spectra in the 15-150 keV band 
were extracted using the tool batbinevt. We applied all the re- 
quired coiTections: we updated them through batupdatephakw 
and generated the detector response matrices using batdrmgen. 
Then we used batphasyserr in order to account for the BAT 
systematics as a function of energy. Finally we grouped the en- 
ergy channels of the spectra by imposing a 3cr (or 2-cr when the 
S/N was too low) threshold on each grouped channel. We fit- 
ted the resulting photon spectra, 0(£') (ph cm^^s^'keV"'), with 
a power law with pegged normalisation (pegpwrlw model un- 
der xspec V. 11.3.2). We extracted several spectra in different 
time intervals: over Tgo, total, spanning the bunch of the first 
three pulses, the fourth pulse alone and that around the peak, 
determined on a minimum significance criterion. The results 
are reported in Table |2] The time-averaged spectral index is 
Py = 1.65 + 0.51 with a total fluence of 5(15 - 150 keV) = 
(3.6 + 0.8) X 10"^ erg cm"^ and a 0.448-s peak photon flux 
of (1 .0 + 0.2) ph cm~^ s~' , in agreement with previous re- 
sults dMarkwardt et al.ir2008 h. The bottom panel of Fig. [T]shows 
marginal evidence for a soft-to-hard evolution: Py passes from 
2.0 + 0.5 (first three pulses) to 1 .4 ± .5 (fourth pulse). 

Following ISakamoto et al.l d2008l) . a GRB is classified as an 
XRR (XRF) depending on whether the fluence ratio 5(25 - 
50 keV)/5(50 - 100 keV) is lower (greater) than 1.32. The flu- 
ence ratio of XRF 080330, 1.5+jJ-^, places it among the XRFs, 
although still compatible with being an XRR burst. Although 
from BAT data alone we could not measure the peak energy 
of the time-integrated vFy spectru m, we tried to fit with a 
smoothed broken power law model (iBand et al ] ll993h by fix- 

^ |http : //nedwww . ipac . caltech . edu/ forms/ calculator . html | 
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Table 1. Best-fit parameters of the 15-150 keV profile decomposed into four pulses using the model by iNorris et alj ( l2005h . 
Uncertainties are Icr. 



Pulse 




A 






w 




k 


Fluence 




(s) 


(10"** erg cm"^ s"') 


(s) 


(s) 


(s) 






(10"** erg cm--) 


1 


0.5 + 0.2 


4.6 ± 0.6 


0.79 + 0.18 


1.30 + 0.26 


2.09 + 0.29 


0.24 


+ 0.16 


8.6+ 1.0 


2 


4.2 ±0.6 


2.0 + 0.3 


1.04 + 0.64 


4.71 + 1.30 


5.75 + 1.44 


0.64 


+ 0.20 


10.8+ 1.9 


3 


7.5 + 0.2 


2.5 + 0.7 


0.32 + 0.19 


0.84 + 0.43 


1.16 + 0.46 


0.45 


+ 0.32 


2.7 + 0.8 


4 


56.2+1.2 


1.3 + 0.2 


2.52 + 1.26 


10.5 + 2.3 


13.0 + 2.4 


0.61 


+ 0.18 


15.5 + 2.3 




20 40 60 

Time from trigger [s] 



100 



Fig. 1. Top panel: 15-150 keV BAT mask-weighted light curve 
(binning time of 0.512 s). The thick solid line shows the result 
of fitt ing the profile with four pulses modelled with Norris pro- 
files (iNorris et al.ll2005l) . Bottom panel: spectral index jSy as a 
function of time. 



ing the low-energy index to -1 dKaneko et al.ll2006l) . given 
that Y - Py + \ > 2 and is very likely dominated by the high- 
energy index j6b . This way we derived the following constraint: 
E-p < 35 keV, in agreement with the upper limit to its rest- 
fra me (intrinsic) va lue, E-pi - E^il + z) < 88 keV, obtained 

bv lRossi et all (I2008.). 

Recently, I Sakamoto et alj ( l2009h calibrated a method aimed 



at estimating E-p from the Ty as measured with BAT, provided 
that 1.3<r<2.3. In the case of XRF 080330, the confidence in- 
terval on r, 2.65 + 0.51, marginally overlaps with the allowed 
range; however, since /Sp is anti-correlated with F and the lower 
limit on F lies within the usable range, we can derive a n upper 
limit to Ep from the relation of Sakamoto et aH (I2009I) . which 
turns out to be 30 keV, in agreement with our previous value. 
These re sults are fully co nsistent with a previous preliminary 
analysis dMarkwardt et al.ll2008.) . 

We constrained Eiso in the rest-frame 1-10^ keV band us- 
ing the upper limit on E^, of 35 keV. Follow ing the prescriptions 
by Amati et al. (2002) and iGhirlanda eTal . (2004), we found 
iiiso < 2.2 X 10^^ ergs. Combined with E-g^i < 88 keV, this places 
XRF 08 0330 in th e Ep,[-E\so space consistently with the Amati 
relation dAmati et al. 2002; Amati 2006.) . 



3.2. X-ray data 

The XRT data were processed using the heasoft package 
(v.6.4). We ran the task xrtpipeline (v.0.11.6) applying cal- 
ibration and standard filtering and screening criteria. Data from 
77 to 1 34 s were acquired in Windowed Timing ( WT) mode and 



the following in Photon Counting (PC) mode due to the faintness 
of the source. Events with grades 0-2 and 0-12 were selected 
for the two modes, respectively. XRT observations went on up 
to 5.9 X 10^ s, with a total net exposure time of 30.6 ks. The 
XRT analysis was performed in the 0.3-10 keV energy band. 

Source photons were extracted from WT mode data in a rect- 
angular region 40 pixels along the image strip (20 pixel wide) 
centred on the source, whereas the background photons were ex- 
tracted from an equally-sized region with no sources. 

Firstly, we extracted the first orbit PC data from 136 to 331 s, 
where the point spread function (PSF) of the source looked un- 
affected by the spacecraft drifting, and extracted the following 
refined position: RA= ll'' 17" 04!68, Dec. ^ -h30° 37' 24^8 
(J2000), with an error radius of 4.0 arcsec (^ Mao & Guidorzil 
2008). We corrected these data for pile-up by extracting source 
photons from an annular region centred on the above position 
and with inner and outer radii of 4 and 30 pixels (1 pixel = 
2'.'36), respectively. The background was estimated from a three- 
circle region with a total area of 30.3 x 10^ pixel^ away from any 
source present in the field. Finally, we re-extracted the source 
photons over the entire first orbit within a larger circular region 
centred on the same position and with a radius of 40 pixels, to 
compensate for the drifting. The light curve of the full first orbit 
data (PC mode) was then corrected so as to match the previous 
one correctly produced in the 136-331 s sub-interval. 

The resulting 0.3-10 keV light curve is shown in Fig. |3] 
(black empty triangles). It was binned so as to achieve a min- 
imum signal to noise ratio (SNR) of 3. The data taken in follow- 
ing orbits were not enough to provide a significant detection and 
only a 3cr upper limit was obtained. 

The X-ray curve can be fitted with a broken power law, 
with the following parameters: \ - 4.8 + 0.4, fb = 163^jq s, 
Q'x,2 = 0.26 + 0.10 Of^/dof = 66/70). The last upper limit clearly 
requires a further break. We set a lower limit on Qrx,3 by connect- 
ing the end of first orbit data with the late upper limit under the 
assumption that the second break occurred at the beginning of 
the data gap. This turned into a-^ -^ > 1.3 (>3cr confidence). The 
later the second break time, the steeper the final decay. 

We extracted the 0.3-10 keV spectrum in two different time 
intervals: i) "XRT-WT" interval, from 77 to 134 s (WT mode), 
corresponding to the initial steep decay; ii) "plateau" interval, 
from 423 to 1507 s (PC mode) corresponding to the following 
flat decay (or "plateau") phase. Source and background spectra 
were extracted from the same regions as the ones used for the 
light curve for the corresponding time intervals and modes. The 
ancillary response files were generated using the task xrtmkarf. 
Spectral channels were grouped so as to have at least 20 counts 
per bin. Spectral fitting was performed with xspec (v. 1 1.3.2). 

We modelled both spectra with a photoelectrically absorbed 
power law (mod el wabs-zwabs-pow), adop ting the photoelectric 
cross section by iMorrison & McCammorj (I1983i) . The first col- 
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umn density was frozen to the weighted average Galactic value 
along the line o f sight to the GRB, A^*f''^ = 1.23 x lO^" cm-^ 
(iKalberla et al.ll20b5) . while the second rest-frame column den- 
sity, A^H.z, was left free to vary. While during the steep de- 
cay we found no evidence for significant rest-frame absorption, 
with a 90% confidence limit of A^h,- < 1.4 x 10^' cm"^, in 
the plateau spectrum we found only marginal evidence for it, 
Nh,z - ^-^-{'^ ^ 10^' cm"^. The spectral index, j6x, varies from 
1.06+°-^^ to O.SO+U I^: the significance of this change is ~ 2.3 cr. 
The best-fit parameters are reported in Table |2] 

3.3. Near-UV/Visible UVOT data 

The Swift UVOT instrument started observing on 2008 March 
30 at 03:42:19 UT, 63 s after the BAT trigger, with a 9.37-s set- 
tling exposure. Since the detectors are powered up during this 
exposure, the effective exposure time may be less than reported. 
We checked the brightness in this exposure with later exposures, 
to confirm that no correction was needed. The first 99.7-s find- 
ing chart exposure started at 03:42:39 UT in the white filter in 
event mode followed by a 399. 8-s exposure in the V filter, also in 
event mode. Due to the loss of lock by the spacecraft star track- 
ers, the attitude information was incorrect. In order to process 
the data, xselect was used to extract images for short time in- 
tervals. The length of the interval was chosen short enough that 
the drift of the spacecraft was mostly within 7", and at most 
14", but long enough to get a reasonably accurate measurement. 
A source region was placed over the position of the source, mak- 
ing checks for consistency with the position of nearby stars, and 
the magnitudes were determined using the ftool uvotevtlc. In 
most cases, an aperture with a radius of 5" was used, and three 
measurements used a slightly larger aperture. No aperture cor- 
rection was made, since the source shape in those cases was very 
elongated. The measured magnitudes were converte d back to the 
origin al count rates, using the UVOT calibration dPoole et akl 
1200 8). These were sub sequently converted to fluxes using the 
method of Poole et al.i (2008), but for an incident power-law 
spectrum with yS=0.8 and for a redshift z= 1 .5 1 . 

3.4. NIR/Visible ground-based data 

Robotically triggered observations with the LT began at 181 s 
leading to the automatic ide ntification by the GRB pipeline LT- 
TRAP ( Guidorzi et al.l l2006) of the optical afterglow at the po- 
sition RA = ll^ 17™ 04^.48, Dec. = +30° 37' 23'.'8 (J2000; 
Icr error radius of 0.2 arcsec). This is consistent within 1.6 cr 
with the refined UVOT position. The afterglow was seen during 
the end of the rise with r' - 17.3 ± 0.1 and subsequently de- 
cay jGomboc et al.ll2008ah . Following two initial sequences of 
3 X 10 s each in the r' filter during the detection mode (DM), 
the multi-colour imaging observing mode (MCIM) in the g'r'i' 
filters was automatically selected. Observations carried on up to 
4.9 ks. 

The FTN observations of XRF 080330 were carried out from 
8.4 to 9.1 hr and again from 31.8 to 33.9 hr with d eep r' and i' fil- 
ter ex posures as part of the RoboNet 1.0 projecfl jGomboc et al.l 
12006 ). 

Calibration was performed against fi ve non-saturated field 
stars with preburst SDSS photometry jCool et al.1 l2008h . by 
adopting their PSF to adjust the zero point of the single images. 
Photometry was carried out using the Starlink GAIA software. 



Magnitudes were con verted into flux densities (mJy) following 
Fukugit a et al.l (Il996h . Results are reported in Table|5] 

Optical iJ-band observations of the afterglow of 
GRB 080330 were carried out with the REM telescope equipped 
with the ROSS optical spectrograph/ima ger on 2008 March 30 , 
starting about 55 seconds after the burst ( D'Avanzo et al.ll2008l) . 
We collected 38 images with typical exposures times of 30, 60 
and 120 s, covering a time interval of about 0.5 hours. Image 
reduction was carried out by following the standard procedures: 
subtraction of an averaged bias frame, division by a normalised 
flat frame. The astrometry was fitted using the USNOBI.CQ 
catalogue. We grouped our images into 18 bins in order to 
increase the signal-to-noise ratio (SNR) and p erformed aperture 
photometry with the SExtractor package dBertin & ArnoutsI 
1 1996) for all the objects in the field. In order to minimise the 
systematics, we performed differential photometry with respect 
to a selection of local isolated and non-saturated standard stars. 

The calibration of NOT images taken with the R filter was 
performed with respect to the converted magnitudes in the R- 
band of the selected set of stars used for the calibration of LT 
and GROND images in the SDSS passbands. We transformed 
the r' and /' magnitudes of the calibration stars dCool et al.l2008l) 
into R and / ma gnitude following the filter transformations of 
iKrisciunas et al.l 11998). 

Hereafter, the magnitudes shown are not corrected for 
Galactic extinction, whilst fluxes as well as all the best-fit mod- 
els are. When the models are plotted together with magnitudes, 
the coiTection for Galactic extinction is removed from the mod- 
els. 

3.5. Spectroscopy 

Starting at a!46 min we obtained a 1800 s spectrum with a low 
resolution grism and a 1.3 arcsec wide slit covering the spectral 
range from about 3500 to 9000 A at a resolution of 14 A with 
the NOT (Fig. 121). The airmass was about 1.8 at the start of the 
observations. The spectrum was reduced using standard meth- 
ods for bias subtraction, flat-fielding and wavelength calibration 
using an Helium-Neon arc spectrum. The rms of the residuals in 
the wavelength calibration were about 0.3 A. The spectrum was 
flux-calibrated using an observation with the same setup of the 
spectrophotometric standard star HD93521. Table |6] reports the 
identified lines. 



4. Multi-wavelength combined analysis 

4.1. Panchromatic light curve 

Figure [3] displays the light curves of the prompt emission (15- 
150 keV) and of the 0.3-10 keV and NIR/visible/UV afterglow 
derived frorn our d ata sets plus some points taken from RAPTOR 
(IWren et al."2008). High -energy fluxes (magnitudes) are refeiTed 
to the right-hand (left-hand) y-axis. First of all, we note that the 
peak time of the last y-ray pulse (Table [U is contemporaneous 
with t he optical fl ash detected by RAPTOR, reported at 58.9 + 
2.5 s (IWren et al.l l2008). 

The initial steep decay observed by XRT is a smooth contin- 
uation of the last y-ray pulse and is thus the tail of the prompt 
GRB emission, and likely to correspond to its high-latitude emis- 
sion. Most notably, during the X-ray steep decay the optical flux 
is seen to rise up to ~300 s and finally a simultaneous plateau is 



^ |http : / /www ■ astro . liv jm . ac . uk/RoboNet/j 



* |http : / /www ■ nof s . navy . mil/data/ f chpix/j 
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Table 2. Best-fit parameters of the 15-150 keV and of the 0.3-10 keV spectra of the y-ray prompt and X-ray afterglow measured 
with BAT and XRT, respectively. The model is an absorbed power law (xspec model: pegpwrlw for BAT and wabs pow for XRT 
data, respectively). Frozen values are reported among square brackets. Also SED modelling results are reported, both derived on the 
broad band from optical to X-rays and optical points alone. The extinction is modelled with a SMC profile as parametrised bv lPeil 
(119921) . 



Interval 


Energy band 
(keV) 


Start time 
(s) 


Stop time 
(s) 




A'h.z 
(lO^' cm-2) 


Mean flux 
(erg cm"^ s"') 




xVdof 


7^90 

ioiai 

Pulses 1-3 
Pulse 4 
Peak 


15-150 

15-150 
15-150 
15-150 


0.0 
— Z.U 
-2.0 
52.9 
0.384 


67.0 
on n 

20.0 
90.0 
0.832 


1.44 + 0.46 

1 -1- n ^ 1 
1 .Oj ± U.J 1 

2.0 + 0.5 
1.4 + 0.5 
1.1+0.6 




(3.3 + 0.8) X 10-^ 
{j.o + V.o) X lU 
(2.2 ± 0.5) X lO-'' 
(1.4 + 0.5) X 10-^ 
(1.0 + 0.2)* 




1.44/6 

J.D/ / 

1.82/7 
5.8/5 
6.1/6 


XRT-WT 
Plateau 


0.3-10 
0.3-10 


77 
423 


134 
1507 


i.06!^:i'^ 
o.80!»:]^ 


< 1.4 


(4.1 +0.3) X 10-"' 
(2.3 + 0.3) X 10-" 




28.3/31 
17.2/24 


SED 2 
SED 3 
SED 4 


opt-X 
opt-X 
opt 


186.8 
423 
78117 


269.4 
1507 
93620 


0.74 + 0.03 
0.79 + 0.01 
0.85 + 0.30 


[2.7] 
2.7 + 0.8 




<0.04 
<0.02 
10+°'" 


8.2/6 
32/34 
3.4/4 


SED 2 
SED 3 
SED 4 


opt 
opt 
opt 


186.8 
423 
78117 


269.4 
1507 
93620 


0.61+0.13 
0.74 + 0.05 
1.05+0.06 






[0] 
[0] 
[0] 


3.1/3 
7.9/8 
4.8/5 



" Rest-frame extinction obtained by modelling the SED with an SMC profile as parametrised bv lPeil ( 119921) . 
* Peak photon flux in units of ph cm"- s"'. 




4000 4500 5000 5500 ^ 6000 6500 7000 

Wavelength (A) 



Fig. 2. The NOT spectrum of the afterglow taken at 46 min post burst. The tick marks and Earth symbols show the absorption 
features (Table |6]l and telluric lines, respectively. The dotted line shows the eiTor spectrum. 



reached at both energy bands, lasting up to ~1500 s, when the X- 
ray observations stopped. This strongly suggests that the plateau 
is emission from a region which is physically distinct from that 
responsible for the prompt emission and its tail (the rapid decay 
phase). 

As discussed in Sect. 14.21 the afterglow does not show evi- 
dence for spectral evolution throughout the observations, except 
for late epochs (t ~ 10^ s), when there is evidence for redden- 
ing. The achromatic nature of the afterglow light curve allows 
for a multi-wavelength simultaneous fit of twelve light curves. 



where only the normalisations are left free to vary independently 
from each other. We consider all of the available passbands: K, 
H, J, z', i', r', V, g\ B, U, UWVl and X-ray, respectively. The 
latter curve is fitted from 300 s onward, so as to exclude the 
initial steep decay. Hereafter we present two alternative combi- 
nations of models, both providing a reasonable description of 
the flux temporal evolution. In both cases we had to add a 2% 
systematics to all of the measured uncertainties to account for 
some residual variability with respect to the models, in order to 
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Fig. 3. Panchromatic light curve of the prompt (BAT, 15-150 keV), X-ray (XRT, 0.3-10 keV) and NIR/visible/UV (FTN, GROND, 
LT, NOT, REM, TAROT, UVOT-W) afterglow o f XRF 0803 30. Magnitudes (high-energy fluxes) are refeiTed to the left-hand (right- 
hand) y-axis. RAPTOR points were taken from I Wren et alj (12008) and shifted by 0.3 mag in order to match contemporaneous r' 
points. The shaded bands indicate the four intervals where we computed the SEDs. 



have acceptable x values and correspondingly acceptable pa- 
rameters' uncertainties. 



4.1 .1 . Multiple smoothly broken power law 

A possible description of the light curves is offered by a multiple 
broken power law (Fig. |4]l. This has the advantage of a more 
straightforward interpretation in terms of the standard fireball 
evolution model due t o synchrotron emi s sion. We started from 
the parametrisation by iBeuermann et al.l (11999!) and added two 
more breaks to finally provide a sufficiently detailed description. 
The fitting function is given by eq. ([1]). 



Fit) = 



(1) 



The free parameters are the normalisation constant (different for 
each curve), Fq, three break time constants, tu (/ = 1,2,3), 
four power-law indices, < a2 < < 04, the smooth- 
ness n. Apart from the normalisations, all the curves share the 
same parameters. Overall, the free parameters and the degrees 
of freedom (dof) total 20 and 184, respectively. Equation ([U 
looks like a piecewise power law only in the following regime: 
tbi ti,2 th3, where each individual term takes over at 
well separate epochs. The light curve of XRF 080330 fits in this 
case, as proven by the best-fit results (first line of Table [3]l and 
shown in Fig. H] The effective break times, f;,i,eff, fi2,eif, fi3,eff, 
i.e. the times at which the model ([T]) changes the power-law 
regime, are simply given by tbien = hi, fi2,eff = (^'j/CP'^^"' 

The goodness of the fit in terms of;t'^/dof is 212/184, cor- 
responding to a non-rejectable P-value of 7.7%. The normali- 
sation constants for the different bands are the following (//Jy): 



Fk = 1077+54, Fh = 942+45, ^ 759+40^ j, ^ 543+26^ ^. ^ 

554+1^, Fr = 464+1^, Fy = 418+^^^^ Fg = 362+}^, Fb = 362 + 68, 
Fu = 277 + 45, Fifvwi = 129 + 52, while the X-ray normal- 
isation, expressed in flux units in the 0.3-10 keV band instead 
of flux density, is F^ = (3.7 + 0.3) x 10"" erg cm^^ ^-i j^e 
effective break times are found be fii,eff - 317 s, fe.eif - 1850 s 
and ?fe3,efF = 82.4 ks, respectively. 

The bottom panel of Fig.|4]shows the residuals of the r' curve 
with respect to the model; the displayed uncertainties do not in- 
clude the 2% systematics added by the fitting procedure. We note 
that between 6 and 7 x 10^ s the model overpredicts the flux by 
2-4cr with respect to the measured values, corresponding to a 
~ 0.1 magnitude difference. However, the later points seem to 
rule out a steeper decay than the modelled one. Alternatively, one 
might interpret this as suggestive of a steeper decay followed by 
a second component thus causing a late flux enhancement. This 
possibility motivated us to provide an alternative description, de- 
scribed in the next Sect. 14. 1.2] 

4.1.2. A two-component model: late-time brightening 

In Figure |5] we modelled the first part (f < lO'* s) with a simple 
smoothly broken power law with a single break time, th\, and two 
power-law indices, ai {02), taking over for at f «c tb\ (t » tb\). 
In order to model the later data points, we had to add a further 
component. A Lorentzian proved successful in this respect, so 
that the complete model used is given by the following equation. 

F{)y Fty 



Fit) 



[(f/fM)""' + itltbiT"']"" 1 + [2 (f - fc) /fw]' 



(2) 



This was used to fit the r' curve. The free parameters are the 
normalisation constant, Fq,,-, the break time, tb\, two power-law 



8 



C. Guidorzi et al.: Prompt and afterglow study of XRF 080330 



10' 



10" 



10' 



^ 10 



10' 



-2 



10' 



10' 



< 



2 

-2 
-4 




5xK 
4xH 
3xJ 
2xz 
1.5xi 
r/R 
0.5xV 
0.3xg 
0.2xB 
O.lxU 
O.OSxUVWl 
3xX(3 keV) 



-H+H : : 1 i 1 1 1 1 H+H : H 

— , — i 


r 1 ^ — > — ' — ' — ^ 

- — ii-^ J : 


_ I ^ ¥ ^ 
" . . 1 1 * . . 1 


1 1 ■ ■ ■ ~ 



10^ 



10^ 10* 10^ 

Observed time from trigger [s] 



10° 



Fig. 4. Top panel: NIR/visible/UV/X-ray light curves expressed in flux densities units, after correction for Galactic extinction. The 
effective wavelength decreases from top to bottom, from K filter all the way down to X-rays. Re d empty circles are public data 
points from Uterature ( Im et al.l l20Q8: Wang et al."2Q08': 'Serg eev et al.l l2008: Mos kvitin et al.lf2008h . Upper Umits are 3cr. The solid 
lines show the same model of a multiple smoothly broken power-law obtained by fitting all the curves simultaneously, allowing for 
different normalisations but with the same fit parameters. Bottom panel: residuals of the r' curve with respect to the corresponding 
model. 



indices, ai and a2, the smoothness n, the Lorentzian normalisa- 
tion, /^L.r, the peak time, t^ and its width, The time-integrated 
flux density of the latter component is n FL,r fw/2. The two terms 
of eq. ^ peak at f,,i = fii(-Qfi/a;2)'''"^"'^'"*^ and tp2 = fc, re- 
spectively. Each of the light curves of the remaining filters were 
fitted with a free scaling factor with respect to the r' curve as 
modelled by eq. (|2]i. The free parameters and the dof total 19 
and 185, respectively. The best-fitresult is shown in Fig.|5] while 
the second line of Table [3] reports the corresponding best fit val- 
ues. The fit is good: ;t'^/dof = 187/185. The scaling factors 
for the remaining bands are the following: fK - 2.31 + 0.12, 
fu = 2.03 + 0.10, /y = 1.65 + 0.09, /, =' 1.40 ± 0.04, = 
1.20 + 0.03, fv = 0.90 + 0.08, fg = 0.78 + 0.02, fe = 0.77^° ;|^, 
fu = 0.59+;j;[^, fuvw\ = 0.28+j;;[j, while the X-ray normaU- 
sation is stifl = (3.7 + 0.3) x 10"" erg cm"^ s"'. The two 
components peak at tpi - 600 s and tp2 - 34.4 ks, respectively. 



We also tried to model the second component with a rising 
and falling smoothly broken power law instead of a Lorentzian. 
However, this brings in too many free parameters, such as the 
slope of the rise, so unless one finds reasons to fix some of them 
to precise values, the fit with such a component turns into highly 
undetermined parameters. 



4.2. Spectral Energy Distribution 

Figure|6]displays four SEDs we derived in as many different time 
intervals (see shaded bands in Fig.|3]l: 

1 . SED 1 includes the last y-ray pulse an d the optical flash de- 
tected by RAPTOR jWren et alj|2008l) . around 60 s; 

2. SED 2 corresponds to the final part of the optical rise, coin- 
ciding with the final part of the X-ray steep decay, spanning 
from 186.8 to 269.4 s; 

3. SED 3 has the broadest wavelength coverage and corre- 
sponds to the plateau phase, from ~400 to ~1500 s. 

4. SED 4 includes NIR/visible measurements around the possi- 
ble late time break in the light curve (Fig.|4]i, at ~10^ s. 

To construct SED 1 w e made use of the RAPTOR measure- 
ment (I Wren et all |2008|) . a UVOT upper limit in the V band 
and of the BAT spectrum of the fourth pulse. Figure |7] displays 
this SED: the solid line shows the best fit with a smoothed bro- 
ken power law used to fit the high-energ y photon spectra of the 
prompt emission of GRBs OBand et al ITl993h . The best-fitting 
parameters are the following: ub - -1.12, jSs = -2.35 and 
Epi = 71 keV (x^/dof = 5.8/5) consistent with the limit on 
Epi derived in Sect. 13.11 The Band indices are photon indices, 
so the corresponding energy indices are 0.12 and 1.35, respec- 
tively. While /3b was constrained by the BAT data themselves, we 
solved the coupled indetermination a^-Ep^i by initially freezing 



C. Guidorzi et al.: Prompt and afterglow study of XRF 080330 



9 




10^ 10^ lO"^ 10^ 10^ 

Observed time from trigger [s] 



Fig. 5. Top panel: same as Fig.|4] In this case the solid lines show the same model of a smoothly broken power-law plus a Lorentzian 
obtained by fitting all the curves simultaneously, allowing for different normalisations but with the same fit parameters. Bottom 
panel: residuals of the r' curve with respect to the corresponding model. 



Table 3. Best-fit parameters of the multi-wavelength fitting procedure of the afterglow light curves. 





(s) 


C2 


tbl 

(s) 




tb3 

(ks) 


Cl'4 


n 


(ks) 


(ks) 


O^Jy) 


A-Vdof 


-0.56™;- 
-0.38!«;i 


3^7+151 
2480!]!'' 


15+'™ 
2 02+''-'*5 


1456- 


1.08 ±0.02 


23.8!^;f 


3 51+"-^" 


•'■^-1.3 

49+0.61 


-iA 4.+ 10.6 


79 7+14.6 
-12.2 


11.9!|^ 


212/184 
187/185 



the low-energy index as to the typical value of - 1 ('K aneko et al.l 
|2006) and then leaving it to vary. The above minimum was so 
found. Although this does not break the degeneracy of both pa- 
rameters (for every Epj < 88 keV there is a value of as for which 
an acceptable fit is given), here it is shown that the extrapolation 
of a typical Band model fitting the spectrum of the last pulse 
matches the optical flux observed during the flash. However, be- 
cause of the lack of measurement of aB from y-ray data, the op- 
tical flux matched by the extrapolation of the Band model may 
still be accidental. 

The time interval of SED 2 corresponds to the first JHK 
GROND frames and spans from 186.8 to 269.4 s (see Fig. [3]l. 
It consists of contemporaneous Vr' JHK frames as well as X- 
rays. The NIR-to-X SED 2 can be fitted with a single power-law 
with j6ox - 0.74+0.03 and negligible dust extinction. The optical 
data alone can be fitted with an unextinguished power law with 
/3o= 0.61+ 0.13 (TableEli. 

SED 3, taken during the plateau, is the richest one includ- 
ing all of the passbands considered in this work, but the y- 
rays (see Fig. [3). Our NIR values are consistent with the JHK 
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Fig. 6. GRB rest-frame SEDs 1 to 4 (shown with asterisks, tri- 
angles, squares and diamonds, respectively). The dashed line 
shows the best-fitting power-law model of SED 2: ySox = 0.79 + 
0.01 and Ay^. < 0.02. X-ray data are not absorption-corrected. 
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Fig. 7. GRB rest-frame SED 1 from observed optical to y-ray 
during the optical flash concomitant with the last y-ray pulse 
at ~ 60 s. The solid line shows the best-fitting smoothed bro- 
ken power law dBand et al.l 1993b with the following parameters: 
cb = -1.1, ySB = -2.35 and E^^i = 71 keV. 
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Fig. 8. GRB rest-frame SED 2 corresponding to the final part of 
the rise of the optical afterglow. 



points of lBloom" &Stani(l2008h . Given the steadiness of the light 
curve and the evidence for no significant colour change along 
the plateau, the SED so obtained is fairly robust. The multi- 
wavelength fitting of the light curves of Figs. l4landl5](Sect. l4. lb 
is dominated by the data points along the plateau phase. Thus, 
we built a SED using the best-fit normalisations (Sects. 1471. H and 
14.1.2b and calculated at the time of 10^ s. We found the same re- 
sults with improved uncertainties, due to the stronger constraints 
imposed by a multi-band fitting. The result of this SED is dis- 
played in Figure |9l The solid line shows the bes t-fit mode l ob- 
tained adopting a rest-frame SMC-extinguished (lPeilll992l) . X- 
ray photoelectrically absorbed power law with jSox - 0.79 + 0.01. 
We note that the point corresponding to the UVWl filter nicely 
agrees with the Lyman absorption at the GRB redshift. The rest- 
frame optical extinction was found to be negligible and a very 
tight limit could be derived. Ay,- < 0.02. Thanks to the more 
precise estimate obtained on /3ox, the estimate of A^h,; improved 
correspondingly: (2.7 + 0.8) x 10^' cm"^. These results are con- 
sistent with what was obtained from the 0.3-10 keV spectrum 
alone and the accuracy of the estimates benefited significantly 
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Fig. 9. GRB rest-frame SED from observed NIR to X-ray during 
the plateau from 400 to 1500 s. The solid line shows the best-fit 
model (SMC profile) corresponding to a single unextinguished 
power law with spectral index /3 - 0.79 ± 0.01 and negligible 
local-frame extinction, Av,z < 0.02 (see text). 



from the inclusion of NIR/visible data. Fitting the optical data 
alone, the result is similar: no need for a significant amount of 
extinction and, more importantly, the same index: /3o = Pox = Px 
(Table |2|, thus ruling out the possibility of a significant redden- 
ing due to dust along the line of sight within the host galaxy. 
As a consequence, two properties are inferred: 

- a negligible dust column density in the circumburst environ- 
ment and along the line of sight to the GRB through the host 
galaxy; 

- a single power-law component accounting for the (observed) 
NIR to X-ray radiation, pointing to a single emission mech- 
anism with no breaks in between. Furthermore, a single 
power-law spectrum implies an achromatic evolution, con- 
sistently with the observations, while the other way around 
is not true. 

The epoch of SED 4 is ~80 ks, i.e. around the final break 
(aftermath of the late brightening) following the light curve de- 
scription given in Sect. 14. 1 .Tl (Sect. 14.1.2b and shown in Fig. H] 
(Fig. |5]l. This includes optical data and is shown in Fig. |6] (di- 
amonds). Data can be fitted either with a single unextinguished 
power-law with j6o - 1.05+0.06 (A v/, fixed to 0) or, alternatively, 
with j6o = 0.85 + 0.30 and some extinction. Ay,, = O.lO+g g^. An 
increase of the extinction along with time seems hard to explain 
physically, so we are led to favour a true reddening at this time. 
Compared with the previous SEDs (Table|2l), SED 4 is redder by 
AjSo - 0.26 + 0.06 (significance of ~6 x 10"^). We point out that 
the reddening is independent of the fit choice, as demonstrated 
from the comparison of the bare power-law indices with no dust 
correction between the earlier and later spectra (Table|2]l. 



5. Discussion 

The 15-150 keV fluence and peak flux of XRF 080330 ai-e typ- 
ical of other XRFs detected by Swift. The X-ray afterglow flux 
places XRF 080330 in the low end of the distribution of the 
GRB s observed by Swift, similarly to the majority of XRFs 
(Sa kamoto et al. I l200i . Moreover, the observed X-ray flux of 
XRF 080330 lies in the low end of both the XRFs sample of 
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Observed time from trigger [s] 

Fig. 10. 0.3-10 keV steep decay curve. The solid line shows the 
best-fitting model of the form F(t)+k (t-toY"^ , where the param- 
eters of F{t) (eq. [U were frozen to their best-fit values obtained 
by the multi-band fit of the afterglow component (Sect. 14.1.1b 
and with k, to and free parameters. Both terms, F(t) and 
the power-law, are shown separately by the dotted and dashed 
curves, respectively. 



Swift considered by ISakamoto et alj (l2008h an d of the XRFs 
sample of BeppoSAX of ID' Alessio et all feOOe). The X-ray af- 
terglow of XRF 080330 does require a remarkable steepening 
after the shallow phase, with a^.s > 13 regardless of the light 
curve modelling (Fig.[3]l. This decay is typical of classical GRBs 
(or steeper), but is in contrast to the fairly shallow decays found 
by ISakamoto et all (12008) for their sample of Swift XRFs. The 
optical flux of the XRF 080330 afterglow is within Icr of the 
distrib ution of the BeppoSAX XRFs sample of ID' Alessio et alj 
(12006') at 40 ks post burst. 

The coincidence of the steep decay observed in the X-ray 
light curve, that is a smooth continuation of the last y-ray 
pulse, suggests that this corresponds to its high-latitude emis- 
sion or the so-called "curvature effect" dFenimore et aTl Il996t 
iKum ar & Panaitescu 2000: ' Dermeij |2004V In the case of a thin 
shell emitting for a short time, the closure relation expected be- 
tween temporal and spectral indices is a - /3 + 2, with the time 
origin to reset to the ejection time of the related pulse, earlier 
than the onset by about 3-4 times the width of the pulse. This 
still holds even if t he emission occurs over a finite range of radii 
dGenet & Granotl2008 ). though in that case the ratio of the ejec- 
tion to pulse onset time difference and the pulse width becomes 
smaller (~ 1 for AT? ~ R). 

We fitted the X-ray decay up to 210 s with the form F(t) + 
k(f - fo) where the parameters of F(t) (eq. [1]) were frozen 
to their coiTesponding best-fit values obtained in Sect. l4.1.Tl and 
the power-law parameters fo, and its normalisation k were 



left free to vary. We obtained: to - 53;^'j'g s and = 2.4;^q5 s 
Cf^/dof = 22/27), as shown by the solid line in Fig.[TOl During 
the steep decay, it is j6x = ^-O^^qq^ (Table |2|. The curvature 
relation is fully satisfied and we note that to does coiTespond 
to the time of the last pulse. Replacing eq. ([T) with eq. (|2]i for 
the underlying component, F(t), the best-fit parameters do not 
change to a noticeable degree. 

The optical afterglow of XRF 080330 exhibited a slow rise 
up to ~ 300 s, followed by plateau out to ~ 2 x lO"' s, after 
which it decayed within a typical power-law index of about 1 . 1 
approximately out to a few 10** s. Then, a sharp break to a decay 
index of 3.5^g3 occuiTed concurrently with an optical redden- 



ing (Sect. 14.1. It Fig. |4]i. Alternatively, after the plateau a more 
gradual transition to a power-law decay index of 2.0 ± 0.8 set in, 
followed by a smooth, red bump (Sect. l4.l"!2l Fig.|5]). We discuss 
each phase separately in the following subsections. 



5.1. Optical afterglow rise 

In the context of the fireball model (e.g. Meszaros 2006 and ref- 
erences therein), the possibility that the peak of the optical af- 
terglow emission coiTesponds to the passage of the peak syn- 
chrotron frequency is ruled out by the lack of spectral evolution: 
y6o should evolve from negative (-1/3) to positive values, while 
we find no evidence for /3o changing before ~8 x 10'* s. 

Another possible interpretation of the optical peak is the on- 
set of the afterglow, a s for ORB 060418 and GRB 060607A 
(Molinari.etiD 120071: 'Jin & Fan' '2007') and, possibly, for 
XRF 071010A as wefl (Covinoetal. 2008). In the case of 
XRF 080330, the rise during the pre-deceleration of the fireball 
within an ISM is much shallo wer than a ~ -3, expected at fre- 
quencies betwee n Vm and Vc (ISari &Piranlll999H G rano tl l2005t 
iJin & Fanll200 7). A wind environment would fit in a better way 
the slow rise of XRF 080330. Under these assumptions and in 
the thin shell case as the duration of the GRB is much shorter 
than the deceleration time, we can estimate the initial Lorentz 
factor, Fo (approximately, twice as large as the Lorentz factor 
at the peak), in a wind-shaped density profile, «(r) = A r^^ (A 
is constant), with s - 2, from the peak t ime and the y-ray ra - 
diated energy, Ei,a dChevalier & Lil(2000t iMoUnari et al.ll2007h. 
For consistency, only the two-component model (Sect. 14.1.21 
Fig. [Sj must be considered. In this case, we take the peak time 
of the first component, tpi = 600 s: assuming 77 = 0.2 (radia- 
tive efliciency), A = 3 x lO-'^ cm"', it turns out Fq < 80, and 
a corresponding deceleration radius smaller than 7 x 10'^ cm. 
As in the case of XRF 071010A, the initial Lorentz factor is 
smaller than those found for classical GRBs. There is no evi- 
dence for the presence of a reverse shock; should the injection 
frequency of the reverse shock lie within the optical passbands, 
it w ould dominate the optical fl ux and exhibit a fast (~r^ ') de- 
cay ('Kobavas hi & Zhangl2007h . not observed here. Nonetheless, 
this can still be t he case if the irijection frequency li es below 
the optical bands (iJin & Fanl2007l: iMundeU et al.l2007h . A weak 
point of this interpretation is that the case s = 2 is ruled out: 
a = s(p + 5)/4 - 3 = 0.79 + 0.01 ( lGranotll2005h . Inverting this 
relation, a value of = 1.4 + 0.1 is required to explain the ob- 
served a = -0.4 + 0.2. This argument, together with the absence 
of reverse shock, whose Fv,max should be much larger compared 
with the forward shock by a factor of ~ F (although see above), 
makes the interpretation of a deceleration through a wind envi- 
ronment somewhat contrived. 

In the context of a single jet viewed off-axis d Granot et alJ 

I2OO5I) . the rising part of the XRF 080330 curve is explained 
by the emission coming from the edge of the jet: as the bulk 
Lorentz factor F decreases, the beaming cone gets progressively 
wider, thus resulting in a rising flux. The peak in the light curve 
is reached when it is F ~ 1 / (0obs - ^o), where 0obs and 6*0 are the 
viewing and jet opening angles, respectively. A ccording to the 
optica l afterglow classification given by IPanaitescu & Vestrandl 
(12008!), XRF 080330 belongs to the class of slow-rising and 
peaking after 100 s events. Those authors found a possible anti- 
correlation between the peak flux and the peak time for a num- 
ber of fast-rising afterglows, followed also by the slow-rising 
class and, in this respect, XRF 080330 is no exception. The sug- 
gested interpretation of the rise is either the pre-deceleration syn- 
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Observed time since trigger [s] 

Fig. 11. r-band afterglow of XRF 080330 compared with the 
XRR 030418 (empty circles), fo r which there is only an up- 
per Umit to its redshift, z < 5 dPuUighan et alJ l2003l). Dat a 
of XRR 030418 have been taken from iRvkoff et all (|2004: 
iFerrero et all (|2003[) : iDulhghan et"aLl (|2003|) . The "so lid line is 
the best fit of the XRF 080330 r curve of Sect. |4Xn while the 
dashed line is the best fit obtained with the same model applied 
to XRR 030418. 



chrotron emission or the emergence of a highly collimated out- 
flow seen ofF-axis. In the latter case, assuming a power-law an- 
gular distribution of the kinetic energy, &{6) oc (9/9^^ (q > 0), 
high values for q correspond to slower rises and dimmer peak 
fluxes, for a fixed off-axis viewing angle (0obs = 26*0). In the 
former case, the anti-correlation is ascribed to different circum- 
burst environment densities for different events: XRF 080330, 
because of the negligible dust extinction, would lie in the high- 
peak flux region, which is not the case. This favours the inter- 
pretation of an ofF-axis jet whose angular distribution of energy 
quickly drops away from the jet axis. 

An example of another XRF whose optical counterpart 
showed a very similar behaviour is XRR 030418. The rise of 
this XRR, for which only an ti pper limit to its redshift (z < 5) 
was obtained (iDullighan et al.ll2003 ), has been interpreted as 
due to the decreasing extinction along with time, caused by 
the dust col umn dens ity crossed by the fireball during its ex- 
pansion (Rv koff et al.l b004). Figure [TT] shows the light curve 
compared with that of XRF 080330. The solid line shows the 
best fit to the r' curve of XRF 080330 of Sect. |4Xn while the 
dashed line shows the same model applied to the XRR 030418 
data. XRR 030418 displays a steeper rise (ai - -1.5), which 
strongly depends on the zero time and could be the same as that 
of XRF 080330 if the time origin was moved by (190 + 50) s 
forward in time (lab frame). However, there is nothing around 
this time in the y-ray light curve of XRR 030418. Apart from 
the different slopes of the rise and the lack of a late-time steep- 
ening in the case of XRR 030418, the plateau and post-plateau 
decay look very similar. If both XRFs are caused by the same 
process, from the spectral (lack of) evolution XRF 080330 dur- 
ing the rise-plateau-initial decay phases we can rule out the de- 
creasing dust column density hypothesis. 

5.2. Plateau 

From the SED extracted around the plateau no break is found 
between optical and X-ray frequencies, with /Sox = 0.79 + 0.01. 
In the regime of slow cooling it is reasonable to assume that 



both optical (vo) and X-ray (v^) frequencies lie between the in- 
jection (Vni) and the cooling (v^) frequencies: v^, < Vq < Vx < 
(Sari et al. 1998). The power-law index of the electron energy 
distribution, p, is given by ySox - ip - yielding p = 

2.58+0.02, fully within the range of values found for other bursts 
(e.g. Starling et al. 2008). The temporal decay index depends 
on the density profile: the ISM (wind) case predicts a value of 
a^3(p- l)/4 = 1.18 + 0.02 (a = 3p/4 - 1/4 = 1.68 + 0.02). 
After the plateau, depending on the light curve modelling, the 
measured decay index is either a = = 1.08+0.02 (Sect. 14.1.11 
Fig.© ora = a2 = 2.0 + 0.8 (Sect. l4T2l Fig.©. W hile the mul- 
tiple smoothly broken power-law description (Sect. 14.1. Il l defi- 
nitely rules out the wind environment, both environments are 
still possible in the two-component model (Sect. 14.1.2b . mainly 
because of the poorly measured decay index , g2- If on e interprets 
the flat decay as due to energy injection dNousek e t al. 200^ 
IZhang et"aLll2006l) . the coiTesponding index would he q ~ 0.3. 

In the ofF-axis jet interpretation, even if we consider an ini- 
tially uniform sharp-edged jet, the shocked external medium at 
the sides of the jet has a significantly smaller Lorentz factor than 
near the head of the jet, and therefore its emission is not strongly 
beamed away from off-beam lines of sight. As a result, either 
an early very shallow rise or decay is expected for a realistic jet 
structure and dynamics (Eichler & Granot 2006). In the case of 
XRF 030723. iGranot et al.l 42005) showed that, for 0obs ~ 2 9o, 
an initial plateau is expected in the light curve. 

Our observations of the afterglow rise of XRF 080 330 rule 
out the interpretation proposed by Yamazaki| (l2009l) of the 
plateau as due to an artifact of the choice of the reference time, 
as all the other rising curves do. 

5.3. Jet break 

According to the light curve description of Sect. 14.1.1! shown 
in Fig. m for which only an ISM environment is possible 
(Sect. 15.2b . after the plateau phase the light curve is expected 
to approach the on-axis light curve with a - 3(p - l)/4. The 
late-time steepening observed around 8 x lO'* s, estimated to be 
Aa = Q'4 - Qf3 -2.4 + 0.4, cannot be produced by the passage of 
the cooling frequency through the optical, as that is expected 
to be as small as Aa - 1/4 (ISM/wind). 

In the off-axis jet interpretation, assuming a value for Oq of a 
few degrees, another advantage of this interpretation is the steep 
late time decay (at > 1 day) as a consequence of joining the post 
jet break on-axis light curve. According to the light curve mod- 
elling given in Sect. l4.1.2l shown in Fig.|5] the gradual steepening 
following the plateau corresponds to the post-jet break emission: 
the observed power-law decay index, 0-2 - 2.0 + 0.8 is consistent 
with the expected a - p - 2.6. We note that the relatively sharp 
jet break favours the ISM environment. Overall, in the context 
of an ofF-axis viewing angle interpretation, the light curve sug- 
gests that 0obs ~ (1.5 - 2)0() as well as an early jet break (at 
< 1 day), which in turn implies a narrow jet with a half-opening 
angle of the order of a few degrees, 0o ~ 0.05. As a simple 
feasibility check, we note that for ^obs < 200 the ratio of the on- 
axis to off-axis Ey,\so is equal to 5^ (assuming that the observed 
energy range includes E-p where most of the energy is radiated), 
where 6 is the ratio of their correspondiiig Dop pler factors and 
therefore of their E-p jEichler & Levinsonl2004l) . In our case, for 
an observed ofF-axis (1 + z}Ep ~ 60 keV, an on-axis value of 
~ 1 MeV would require 6^1 + [ro(6'obs - fo)]^ ~ 17 and 
ro(0obs - Oo) ~ 4, which for flobs - ~ (0.5 - l)9o and 9o ~ 0.05 
gives To ~ 80 - 160. Here Fq is the initial Lorentz factor at the 
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Table 4. Summary of the main properties of XRF 080330. 



Name 


Value 


z 


1.51 


5(15 - 150 keV) 


(3.6 ± 0.8) X 10"' erg cm"- 


P(15 - 150 keV) 


(1.0 ± 0.2) ph cm"- s"' 




< 35 keV 


£p,i =£p(l+z) 


< 88 keV 


£iso (15-150 keV, obs frame) 


(2.1 ±0.5) X 10^' ergs 


fi'iso (1-10^ keV, GRB frame) 


< 2.2 X 10'^ ergs 


fje( (obs frame) 


< 1 day 


6q (jet opening angle)'"' 


few degrees 


fobs (viewing angle)*"* 


-(1.5-2)00 



Under the assumption of a single off-axis jet. 



edge of the jet. More realistically, the jet would not be perfectly 
uniform with extremely sharp edges, and instead Fq is expected 
to be lower at the outer edge of the jet and larger near its center 
(where it could easily reach several hundreds in our illustrative 
example here). In this case 6^ ~ 300 so that the observed Ziy^so 
in the 15 - 150 keV range, which is 2 x 10^' erg, would imply 
an on-axis value for Ey^iso of ~ 10^^ erg, which for a narrow jet 
with 0obs ~ 0.05 would correspond to a true energy of the order 
of 10^' erg. This demonstrates that this scenario can work for 
reasonable values of the physical parameters. We point out that 
the estimate of the on-axis iiy iso of ~ lO^'* erg is for a wide en- 
ergy range containing Ep, since in our illustrative example most 
of the energy is released within the observed range. A more ac- 
curate estimate of the break time and of the corresponding open- 
ing angle is difficult, due to the degeneracy involved in the light 
curve modelling (Figs. |4] and |5]l. Table |4] summarises the main 
properties of XRF 080330. 

5.4. Late time red bump 

Overall, the two-component description of the Ught curves of 
Sect. 14. 1 .2] shown in Fig. |5] appears to be slightly favoured over 
the multiple smoothly broken power-law of Fig.|4] So, irrespec- 
tive of the nature of the rise and plateau, we speculate on the pos- 
sible nature of the second component, modelled in Fig.|5]as a late 
time bump. Clearly, a SN bump, su ch as that possibly observed 
in the light curve of XRF 030723 dFvnbo et al.ll2004l) . is ruled 
out mainly because that is expected to peak days later, which is 
incompatible with one single day after XRF 080330; not to men- 
tion the too high redshift of XRF 080330 for a 1998bw-Hke SN 
to be detected. 

Alternatively, a density bump seems a viable solution, given 
that Vo < Vc (e.g. Lazzati et al. 2002; Guidorzi et al. 2005), 
although the explanation of the observed contemporaneous 
reddening requires ad hoc assumptions, such as the case of 
GRB 050721, which showed similar properties to XRF 080330 
(same jSox with no breaks between optical and X-ray, late time 
redder optical bump). In that case, the observed reddening was 
explained as due to the presence of very dense clumps sur- 
viving the GRB radiation and with a small covering factor 
jAntonelli et al.ll2006,) . 

If the late reddening is due to the passage of through 
the optical bands, in addition to what argued in Sect. 15.31 an- 
other weak point of the multiple smoothly broken power-law 
description of Fig. |4] (Sect. 14.1.1b is the chromatic change of 
A^o = 0.26 + 0.06 we observe in the optical bands around 



8 X 10 s. The passage of the cooling frequency does not ex- 
plain it: the observed reddening would be 0.5, i.e. twice as much. 
This could still be the case, if our measurement might have taken 
place in the course of the spectral change, as the broken power- 
law spectrum is a simple approximation. However, since it is 
Vc > Vx during the plateau because of the unbroken power-law 
spectrum between optical and X-ray, Vc would have decreased 
very rapidly, thus making this option not reasonable: if at 10^ s 
it is Vc > Vx = 10'^ Hz, at 10^ s it should be Vc > 10" Hz » Vo, 
because Vc oc r'^^ (ISM case), thus this possibility is to be ruled 
out. Likewise, in the wind case it is oc f'^^. The times and 
tg, at which would cross the most redward and blueward fil- 
ters, K and g', would differ by a factor of (v^/v^)^ ~ 20, which 
looks incompatible with the light curves. Furthermore, the ob- 
served reddening rules out the wind case, as /3o should decrease 
from p/2 to (p - l)/2. 

Although an energy injection to the blast wave (forward 
shock) can explain the bump feature in the light curve, it is 
difficult to explain the reddening if we consider only the for- 
ward shock emission, as we have discussed. A possible expla- 
nation for the reddening is that the rebrightening is due to the 
short-lived (At ~ t) reverse shock of a slow shell which caught 
up with the shock front and increased its energy through a re- 
freshed shock (e.g. Kumar & Piran 2000; Granot et al. 2003; 
Johannesson et al. 2006): since that shock is going into a shell of 
ejecta, rather than the external medium, it can have a much larger 
es (magnetised fireball: Zhang et al. 2003; Kumar & Panaitescu 
2003; Gomboc et al. 2008b) and, therefore, a lower Vc, quite 
naturally; for an ISM where Vc decreases with time, then Vc 
of the reverse shock could be around the optical for reason- 
able model parameter values. The need for a separate compo- 
nent to explain a chromatic break in th e light curve was als o 
suggested in the case of GRB 061126 dGomboc et al.ll2008bl) . 
Notably, the final steepening after 10^ s, which in the mod- 
elling of Sect. l4.1.TI (Fig. |4]| is described with - 3.5;^q3, is 
compatible with the high-latitude emission of the reverse shock: 
a - ;0o,iate + 2 - 3.05 +0.06. The jet break might happen slightly 
earlier than the break time in the optical light curve. 

A somewhat more contrived way to explain the bump is 
the appearance of a second narrower jet in the two-componen t 
jet model, as proposed for XRF 030723 (iHuang et al .1120041) . 
In this model, the viewing angle, is within or slightly off the 
cone of the wide jet and outside the narrow jet. The plateau 
phase would reflect the deceleration of the wide jet ( Granot et al.l 
2006). Depending on the isotropic-equivalent kinetic energy of 
the wide and narrow jet, iSw.iso and /Sn.iso, on the jet opening 
angles, 6o^„ and 6o^„, on the initial bulk Lorentz factors, yo.w 
and 7(),n, respectively, as well as on the viewing angle, the af- 
terglow emission of either comp onent is dominant a t different 
times. According to the results of iHuang et al.l (|2004|) . the light 
curve of XRF 080330 could be qualitatively explained as fol- 
lows: the first component obtained in Sect. 14. 1.2] (Fig. 15) repre- 
sents the contribution of the wide jet dominating at early times: 
the rise could be due either to the afterglow onset (in a wind en- 
vironment) or to a viewing angle slightly beyond the wide jet 
opening angle: ^obs Z Gq w. The appearance of the second com- 
ponent would mark the deceleration and lateral expansion of the 
narrow jet, the peak time corresponding to the case yn ~ I /dobs- 
Unlike XRF 030723, which showed a relatively sharp late-time 
peak, the bump exhibited by XRF 080330 looks less sharp and 
pronounced. Although this might suggest a relatively lower en- 
ergy of the narrow jet compared to XRF 030723, yet we cannot 
exclude the case ^w.iso ^ ^n.iso. The ratio of the observed en- 
ergies of XRF 080330, of about 0.6 according to the modelling 
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of Sect. 14.1.21 corresponds to the ratio between the early and 
the late time emissions of the wide and narrow jets, respectively. 
Depending on the values of ^obs, ^o.w, 6'o,n, a comparable ratio of 
observed energies, such as that observed for XRF 080330, can 
still be obtained in the case of a much more energetic narrow jet. 

Such a model turned out to be successful in accounting for 
the naked-eye GRB 0803 19B (Racusin et al...2008) : in that case, 
the two coUimation-corrected energies were comparable, while 
the isotropic-equivalent energy of the narrow jet (0o,n - 0.2°) 
was about 400 times larger than that of the wide jet (0o,w = 4.0°). 

However, in the context of the two-component model the ex- 
planation of the late time reddening simultaneously with the ap- 
pearance of the narrow component emission requires two dif- 
ferent values of p for the two jets, which does not look reason- 
able on a physical ground. Another option is that the cooling fre- 
quency of the second jet might be around the optical at the time 
of the bump: however, if the shock micro-physics parameters {p, 
eg, £(,), are the same for the two shocks, as expected on phys- 
ical grounds, and obviously the external medium is the same, 
then the only thing that is different as far as Vc is concerned, 
is ^iso. Since the presence of the bump requires iin.iso > £^w.iso, 

1 /2 

then this would not work for the wind case, where Vr oc E.. 

^ ISO 

— 1/2 

Even in the ISM case, where Vc oc £. ' the fact that Vc > Vv 

' ^ ISO ' C A 

at 10^ s, and therefore Vc > 10'^ Hz at 10^ s for the wide jet, 
requires ^n.iso/^^w.iso ^10^ (which is very extreme) in order for 
Vc.nClO^ s) to be around 10'^ Hz (required by the observed red- 
dening). Therefore, while this could work in principle, in prac- 
tise it requires extreme parameters. In particular, the amplitude 
of the bump suggests that ~ E^,, and therefore from the re- 
quired -En.iso/^w.iso Z 10"* it would follow On/O^j < 10"^, which 
seems pushed to the extreme. 

Overall, these considerations make the single off-axis jet in- 
terpretation much more compelling. 



6. Conclusions 

XRF 080330 is representative of the XRR and XRF classes of 
soft GRBs. Its y- and X-ray properties of both prompt and high- 
energy afterglow emission place it in the low-flux end of the 
distribution. The multi-band (NIR through UV) optical curve 
showed an initial rise up to ~300 s, followed by a ~2-ks long 
plateau, temporally coinciding with the canonical flat decay of 
X-ray afterglows of all kinds of GRBs, followed by a gradual 
steepening and a possible jet break. We provided two alterna- 
tive descriptions of the light curve: a piecewise power-law with 
three break times, the last of which occurring around 8 x 10"* s 
and followed by a sharp steepening, with the power-law decay 
index changing from 1.1 to 3.5^^ 3. The SED from NIR to X-ray 
wavelength is fitted with a simple power-law with ySox - 0.8 and 
negligible GRB-frame extinction. Ay < 0.02 adopting a SMC- 
like profile, with no evidence for chromatic evolution during the 
rise, plateau and early (< 8 x lO'* s) decay phases. However, af- 
ter the possible late time break we observe a reddening in the 
optical bands of A/3o - 0.26 + 0.06, which cannot be accounted 
for in terms of the synchrotron spectrum evolution of a standard 
afterglow model, unless a different description of the light curve 
is considered. In the alternative model of the light curves, we 
identified two distinct components: the first is modelled with a 
smoothly broken power law and fits the rise plateau and early 
decay of the afterglow, while the second, taking over around 
8 X 10'* s, is modelled with an energy injection episode peak- 



ing at 34^g' ks and with a time-integrated energy of ~6Q% that 
of the first component. 

The X-ray light curve consists of the initial steep decay, 
which is likely the high-latitude emission of the last y-ray pulse. 
At the same time, the optical afterglow rises up to a plateau, tem- 
porally coincident with the X-ray flat decay. In this case, we col- 
lected strong evidence that the emission mechanism during this 
phase is the same from optical to X-rays and is consistent with 
synchrotron emission of a decelerating fireball with an electron 
energy distribution power-law index of p = 2.6. 

The lack of spectral evolution throughout the rise, plateau 
and early decay argue against a temporally decreasing dust col- 
umn density claimed to explain similar optical light curves of 
past soft bursts. 

The optical rise (a ~ -0.4) is too slow for the afterglow onset 
within a uniform circumburst medium, but could still be the case 
if a wind environment is considered. In this case, under standard 
assumptions we constrained the Lorentz factor of the fireball to 
be smaller than 80, thus confirming the scenario of XRFs as less 
relativistic GRBs. However, we found that the interpretation of 
a single-component off-axis jet with an opening angle of a few 
degrees and a viewing angle about twice as large, can explain the 
observations: this not only accounts for the light curve morphol- 
ogy, but also explains the soft nature of XRF of the y-ray prompt 
event. The reddening observed at 8 x 10"* s can be interpreted as 
the short-lived reverse shock of an energy injection caused by 
a slow shell which caught up with the fireball shock front, also 
responsible for the contemporaneous bump in the light curve. A 
two-component jet could also work, but would introduce more 
free parameters and would require extreme conditions. 

The interpretation of the late bump as produced by a den- 
sity enhancement in the medium swept up by the fireball cannot 
be ruled out, although the reddening seems to require ad ho c 
explanations. In this case, as shown bv lNakar & Granoj (l2007h . 
it is hard to produce a flux enhancement with density inhomo- 
geneities, although it is not excluded given the lack of a sharp 
rise in this bump. 

Overall, the XRF 080330 optical and X-ray afterglows 
properties have also been observed in many other GRBs 
(Panaitescu & Vestrand 2008). This both supports the view of a 
common origin of XRFs and classical GRBs, which form a con- 
tinuum and do not call for distinct mechanisms. The importance 
of a prompt multi-wavelength coverage of the early phases of a 
GRB is clearly demonstrated in the case of XRF 080330. 
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Table 5. Optical photometric set of XRF 080330. 



Mid time 


Exposure 


Magnitude^'"' 


Filter 


Telescope 


(s) 


(s) 








438 


30 


17.68 ± 0.03 


SDSS-g 


LT 


695 


30 


17.74 ± 0.04 


SDSS-g 


LT 


963 


60 


17.75 ±0.03 


SDSS-g 


LT 


1323 


60 


17.82 ±0.04 


SDSS-g 


LT 


1677 


60 


17.98 ±0.03 


SDSS-g 


LT 


2036 


60 


18.07 ±0.03 


SDSS-g 


LT 


2385 


60 


18.27 ± 0.05 


SDSS-g 


LT 


2737 


60 


18.36 ±0.04 


SDSS-g 


LT 


3081 


60 


18.50 ± 0.05 


SDSS-g 


LT 


3426 


60 


18.57 ± 0.05 


SDSS-g 


LT 


3832 


120 


18.76 ± 0.03 


SDSS-g 


LT 


186 


10 


17.65 ± 0.05 


SDSS-r 


LT 


207 


10 


17.60 ± 0.03 


SDSS-r 


LT 


228 


10 


17.51 ± 0.03 


SDSS-r 


LT 


295 


10 


17.39 ± 0.05 


SDSS-r 


LT 


317 


10 


17.34 ± 0.04 


SDSS-r 


LT 


339 


10 


17.37 ± 0.03 


SDSS-r 


LT 


609 


30 


17.38 ± 0.04 


SDSS-r 


LT 


864 


30 


17.43 ± 0.04 


SDSS-r 


LT 


1204 


60 


17.44 ± 0.05 


SDSS-r 


LT 


1559 


60 


17.58 ± 0.03 


SDSS-r 


LT 


1918 


60 


17.72 ± 0.04 


SDSS-r 


LT 


2268 


60 


17.86 ± 0.04 


SDSS-r 


LT 


2622 


60 


17.98 ± 0.03 


SDSS-r 


LT 


2967 


60 


18.10 ±0.03 


SDSS-r 


LT 


3312 


60 


18.22 ±0.05 


SDSS-r 


LT 


3659 


60 


18.35 ±0.03 


SDSS-r 


LT 


4579 


3 X 10 


18.65 ±0.04 


SDSS-r 


LT 


4688 


30 


18.65 ±0.04 


SDSS-r 


LT 


30383 


200 


20.78 ± 0.06 


SDSS-r 


FTN 


30605 


200 


20.89 ± 0.06 


SDSS-r 


FTN 


30828 


200 


20.90 ± 0.06 


SDSS-r 


FTN 


31050 


200 


20.91 ± 0.06 


SDSS-r 


FTN 


31271 


200 


20.92 ± 0.06 


SDSS-r 


FTN 


31494 


200 


20.84 ±0.06 


SDSS-r 


FTN 


31717 


200 


20.89 ± 0.06 


SDSS-r 


FTN 


31940 


200 


20.79 ± 0.06 


SDSS-r 


FTN 


32163 


200 


20.86 ± 0.07 


SDSS-r 


FTN 


32385 


200 


20.71 ±0.07 


SDSS-r 


FTN 


116557 


20 X 200 


23.3 ± 0.2 


SDSS-r 


FTN 


523 


30 


17.22 ± 0.04 


SDSS-i 


LT 


779 


30 


17.20 ± 0.05 


SDSS-; 


LT 


1080 


60 


17.26 ± 0.04 


SDSS-; 


LT 


1439 


60 


17.38 ± 0.03 


SDSS-; 


LT 


1799 


60 


17.50 ± 0.04 


SDSS-; 


LT 


2155 


60 


17.64 ± 0.05 


SDSS-; 


LT 


2503 


60 


17.75 ± 0.03 


SDSS-! 


LT 


2850 


60 


17.88 ± 0.04 


SDSS-; 


LT 


3195 


60 


17.98 ± 0.05 


SDSS-; 


LT 


3540 


60 


18.12 ±0.05 


SDSS-! 


LT 


4908 


30 


18.56 ± 0.05 


SDSS-; 


LT 


64063 


4x 150 


21.27 ± 0.09 


SDSS-; 


LT 


120821 


10x200 


> 21.0 


SDSS-; 


FTN 


2890 


1476 


18.24 ± 0.02 


SDSS-g 


GROND 


4699 


1414 


18.88 ±0.02 


SDSS-g 


GROND 


6006 


459 


19.20 ± 0.02 


SDSS-g 


GROND 


6786 


459 


19.40 ± 0.02 


SDSS-g 


GROND 


85843 


11522 


22.60 ± 0.02 


SDSS-g 


GROND 


171989 


5866 


> 24.46 


SDSS-g 


GROND 


254165 


5655 


> 24.69 


SDSS-g 


GROND 


2890 


1476 


18.05 ± 0.02 


SDSS-r 


GROND 


4699 


1414 


18.65 ±0.02 


SDSS-r 


GROND 


6006 


459 


18.99 ±0.02 


SDSS-r 


GROND 


6786 


459 


19.18 ±0.02 


SDSS-r 


GROND 


85843 


11522 


22.29 ±0.02 


SDSS-r 


GROND 


171989 


5866 


> 24.15 


SDSS-r 


GROND 


254165 


5655 


> 24.33 


SDSS-r 


GROND 


2890 


1476 


17.86 ±0.02 


SDSS-! 


GROND 


4699 


1414 


18.43 ±0.02 


SDSS-; 


GROND 


6006 


459 


18.77 + 0.02 


SDSS-i' 


GROND 


6786 


459 


18.95 ±0.02 


SDSS-! 


GROND 


85843 


11522 


22.00 ± 0.03 


SDSS-I 


GROND 
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Table 5. continued. 



Mid Time 


Exposure 


Mag^''' 


Filter 


Telescope 


(s) 










171989 


5866 


> 23.28 


SDSS-i 


GROND 


254165 


5655 


> 23.47 


SDSS-i 


GROND 


2890 


1476 


17.69 ± 0.02 


SDSS-z 


GROND 


4699 


1414 


18.27 ± 0.02 


SDSS-z 


GROND 


6006 


459 


18.58 ± 0.02 


SDSS-z 


GROND 


6786 


459 


18.77 ± 0.03 


SDSS-z 


GROND 


85843 


11522 


21.82 ± 0.04 


SDSS-z 


GROND 


1 T1 non 


5866 


> 22.68 


SDSS-z 


CjKUJND 




JOJJ 


> ZJ. lO 


CT-\CC ^ 

ouoo-z 




ZZO 


oU 


1 OA _i_ A 1 A 

io.yu ± u. iu 




LrKwJMJ 


3iD 


oU 


io.yz ± u. lu 


f(c) 


OKUJNL) 


jjy 


60 


io.ou ± U. lU 


J ' 


OKUJNL) 


/oo 


60 


io.oo ± U. iU 


7(c) 


LjKUJNL) 


871 


60 


16.85 ± O.IU 


lie) 


CjRONU 


972 


60 


16.92 ± O.lU 




CjKONU 


1130 


120 


16.93 ± O.IU 


7(C) 


CjKUNU 


1^1 / 


1 on 

iZU 


1 7 n4 -1- n 1 n 

1 / .U'4 It U. iU 


J 


OKWiMl-^ 


1DU4 


1 OA 

IzU 


1 T AA _i_ A 1 A 

1 /.UU ± U.iU 


J 


OKUiNU 


TOO A 

zyzu 


IzUU 


n/lQj-AIA 

1 /.4o ± U.IU 


7(C) 


OKUJNU 


4/jO 


IzUU 


1 Q AO _i_ A 1 A 

lo.Uz ± U.IU 


7(C) 


LiKUINU 






1 Q ayl j_ A 1 A 
lo.j4 ± U.IU 


7(c) 


OK WIN Lf 


OJOOO 


yyou 


01 j_ A 1 1 
Zl.Dl ± U.l 1 


7(C) 


LiKUiNU 


1 /zui J 


DIOU 


^ 01 QO 

> zi.yy 


7(C) 


LiKUiNU 




4oUU 


00 1 ^ 

> ZZ. 1 J 


7(c) 


OKWiMU 


228 


60 


16.78 ± U.IU 


17(C) 


LjKUNU 


336 


60 


1 C TA 1 A 1 A 

16. /U ± U.IU 


17(c) 


LjKONU 


448 


60 


1 *1 *1 C 1 A 1 A 

16.65 ± U.IU 


77(C) 


LjKUJND 


559 


60 


1 *1 TA 1 A 1 A 

16. /O ± O.IU 


77(C) 


CjKUJNU 


00/ 


60 


1 A AO _i_ A 1 n 

lo.oZ ± U. iU 


77(c) 


LrKUJNU 


/Oo 


60 


1 A A O _i_ ("^ in 

lO.Oo ± U. iU 


77(c) 


LjKUJNU 


o/l 


oU 


1 A A ^ _i_ ("A 1 I'A 

lo.o5 ± U. iU 


77(c) 


LrKUJMJ 


y /z 


OU 


1 A TA _i_ A 1 A 
10. /U ± U.iU 


77(c) 
77' ' 


OKUINU 


1 1 ■2A 


izU 


1 A 'TA _i_ n 1 n 
Id. /U ± U.IU 


77(c) 




1 "2 1 T 


1 OA 

IzU 


1 A TA _i_ A 1 A 

10. /U ± U.IU 


t7(c) 


OKUJNU 


1 JLW- 


1 on 

iZU 


1 A HA 4_ n 1ft 
lO. /*+ ± U. iU 


f7(c) 


VjivWIN lJ 


loyu 


1 OA 

IzU 


1 A Q< _i_ A 1 A 

lO.oD ± U.IU 


t7(c) 


OKUiNlJ 


zyzu 


1 OAA 

IzUU 


1 T Oyl j_ A 1 A 
1 /.Z4 ± U.IU 


7j(c) 


LiKUINU 


4/ JO 


1 OAA 

IzUU 


1 T TQ _i_ A 1 A 

1 /. /o ± U.IU 


7j(c) 


LiKHJINU 


f^fWA 
0Uo4 


Q^A 

yju 


1 Q AA j_ A 1 A 

18. UO ± U. lU 


77(c) 


LrKUJNU 


OJOOo 


yyou 


01 1A_i.A11 
Zl . 10 ± U. 1 1 


77(c) 
77' 


LrKUlN U 


i / zU i J 


J lOU 


> Zl . JO 


77(c) 
72' ' 


LrKUINU 


ZD4zUU 


4oUU 


O 1 OQ 

> Zl .Z8 


77(c) 
77 


LrKUINU 


ZziS 


60 


1 A "71 _i_ n 1 A 

lo. /i ± U. lU 


f (c) 


CjKUJNL) 


3jO 


oU 


10.4 / ± U. lU 


t^(c) 




448 


60 


16.56 ± 0.10 


^(c) 


LiKOJND 


559 


60 


1 A /I O 1 1 /^^ 

16. 4Z ± U.IU 


7'(C) 

A' ' 


LjKUNU 


OO/ 


oU 


1 A _i_ A 1 A 

Id. J / ± U.IU 


7^(C) 

A' ' 


(jKUJNU 


768 


60 


\ C AC 1 A 1 A 

1d.4o ± U.IU 


f(c) 

A*- -* 


CjRuINU 


871 


60 


\ C AC 1 A 1 A 

16.45 ± U.IU 


f(c) 

A*- -* 


CjKuINU 


972 


60 


1 A CC 1 A 1 A 

16.55 ± U.IU 


f(c) 

A*- -* 


CjKUiNU 


131 / 


1 on 

IZU 


1 A ylC _i_ A 1 A 


jt'(c) 

A 


VjlvWlN LJ 


1504 


120 


16.53 ± 0.10 


^(c) 


OROND 


1690 


120 


16.77 ±0.10 




GROND 


2920 


1200 


1 T AO . A l A 

l/.Oo ± 0.10 


f[(c) 


GROND 


4736 


1200 


17.62 ± 0.10 


/f(c) 


GROND 


85868 


9960 


20.93 ±0.13 


/f(c) 


GROND 


172013 


5160 


> 20.37 


/f(c) 


GROND 


254200 


4800 


> 20.49 


KC^) 


GROND 


1238 


60 


17.33 ±0.02 


R 


NOT 


1609 


300 


17.46 ± 0.01 


R 


NOT 


1927 


60 


17.61 ± 0.01 


R 


NOT 


89.4 


30 


17.90 ±0.19 


SDSS-r equivalent 


REM 


167.6 


30 


17.61 ±0.13 


SDSS-r equivalent 


REM 


246.2 


30 


17.41 ±0.12 


SDSS-r equivalent 


REM 


354.2 


30 


17.34 ±0.12 


SDSS-r equivalent 


REM 


432.9 


30 


17.40 ±0.11 


SDSS-r equivalent 


REM 


511.1 


30 


17.42 ±0.13 


SDSS-r equivalent 


REM 


590.5 


30 


17.41 ±0.13 


SDSS-r equivalent 


REM 


663.6 


30 


17.42 ± 0.08 


SDSS-r equivalent 


REM 
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Table 5. continued. 
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Mid Time 


Exposure 


Mag"'" 


Filter 


Telescope 


(s) 










732.7 


30 


17.48 ± 0.09 


SDSS-r equivalent 


REM 


802.7 


30 


17.47 ±0.10 


SDSS-r equivalent 


REM 


871.8 


30 


17.40 ± 0.09 


SDSS-r equivalent 


REM 


971.1 


60 


17.45 ± 0.05 


SDSS-r equivalent 


REM 


1099.9 


60 


17.50 ± 0.06 


SDSS-r equivalent 


REM 


1403.1 


30 


17.65 ±0.11 


SDSS-r equivalent 


REM 


1471.4 


30 


17.57 ±0.10 


SDSS-r equivalent 


REM 


1540.5 


30 


17.62 ±0.11 


SDSS-r equivalent 


REM 


1640.7 


60 


17.58 ± 0.06 


SDSS-r equivalent 


REM 


1769.5 


60 


17.73 ± 0.07 


SDSS-r equivalent 


REM 


50.4 


60.0 


> 16.6 


clear 


TAROT 


120.9 


67.8 


17.62 ± 0.60 


clear 


TAROT 


194.1 


66.6 


17.38 ± 0.46 


clear 


TAROT 


249.3 


30.6 


16.96 ±0.40 


clear 


TAROT 


332.7 


89.4 


17.17 ±0.25 


clear 


TAROT 


429.0 


91.2 


17.07 ± 0.20 


clear 


TAROT 


528.7 


90.0 


16.83 ± 0.29 


R 


TAROT 


626.4 


90.0 


16.83 ± 0.22 


clear 


TAROT 


723.0 


90.0 


16.86 ± 0.22 


clear 


TAROT 


927.0 


90.0 


17.52 ± 0.46 


clear 


TAROT 


1023.9 


89.4 


17.15 ±0.30 


clear 


TAROT 


1123.9 


90.0 


17.00 ± 0.39 


R 


TAROT 


1222.5 


89.4 


16.73 ± 0.21 


clear 


TAROT 


1318.5 


89.4 


16.75 ± 0.25 


clear 


TAROT 


1418.3 


90.0 


> 16.8 


R 


TAROT 


68 


10 


> 17.38 


V 


UVOT 


103 


40 


20.06 ± 0.397 


white 


UVOT 


128 


10 


18.55 ± 0.20 


white 


UVOT 


138 


10 


18.36 ± 0.44 


white 


UVOT 


148 


10 


18.28 ± 0.39 


white 


UVOT 


158 


10 


18.39 ± 0.43 


white 


UVOT 


168 


10 


18.22 ± 0.38 


white 


UVOT 


178 


10 


18.05 ± 0.32 


white 


UVOT 


214 


50 


17.74 ±0.19 


V 


UVOT 


264 


50 


17.72 ± 0.22 


V 


UVOT 


314 


50 


17.75 ± 0.26 


V 


UVOT 


364 


50 


17.52 ± 0.23 


V 


UVOT 


414 


50 


17.37 ±0.17 


V 


UVOT 


464 


50 


17.79 ± 0.24 


V 


UVOT 


514 


50 


17.44 ±0.18 


V 


UVOT 


557 


36 


17.37 ± 0.20 


V 


UVOT 


629 


20 


> 18.02 


UVWl 


UVOT 


654 


20 


17.14 ±0.16 


U 


UVOT 


673 


10 


17.81 ±0.23 


B 


UVOT 


688 


10 


17.76 ±0.14 


white 


UVOT 


733 


20 


17.74 ± 0.34 


V 


UVOT 


782 


20 


> 18.04 


UVWl 


UVOT 


788 


172 


> 18.77 


UVW2 


UVOT 


807 


20 


16.66 ±0.17 


U 


UVOT 


827 


10 


17.66 ± 0.22 


B 


UVOT 


838 


10 


17.55 ±0.12 


white 


UVOT 


848 


10 


17.91 ±0.15 


white 


UVOT 


992 


796 


> 18.55 


UVM2 


UVOT 


1010 


50 


17.37 ±0.17 


V 


UVOT 


1060 


50 


17.37 ±0.17 


V 


UVOT 


1110 


50 


17.41 ±0.17 


V 


UVOT 


1160 


50 


17.72 ± 0.21 


V 


UVOT 


1210 


50 


17.36 ± 0.21 


V 


UVOT 


1260 


50 


17.75 ± 0.26 


V 


UVOT 


1310 


50 


18.11 ±0.29 


V 


UVOT 


1360 


51 


17.77 ± 0.22 


V 


UVOT 


1426 


10 


17.67 ± 0.29 


UVWl 


UVOT 


1475 


20 


18.14 ±0.22 


B 


UVOT 


72100 


5887 


> 20.10 


V 


UVOT 


307644 


47917 


> 21.49 


B 


UVOT 


307940 


47715 


> 20.53 


V 


UVOT 


307583 


47956 


> 21.15 


U 


UVOT 


307483 


48074 


> 21.58 


UVWl 


UVOT 


391039 


40970 


> 23.32 


white 


UVOT 



" Values are not corrected for Galactic extinction. 

Errors at the 68% confidence level and upper limits (3 cr) are given. 
AB magnitudes. 
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Table 6. Lines detected in the XRF 080330 spectrum. 



A^u^ [Al 


"■resi 




Feature 


obs ^ ' 


3889.4 


1548.2/1550.8 


1.5101 


Civ 


4.3±0.7 


4028.8 


1608.5 


1.5047 


Fen 


2.0±0.7 


4198.2 


1670.8 


1.5127 


Ain 


2.0±0.7 


5093.5 


2026.1 


1.5139 


Znii 


4.3+0.4 


5190.3 


2062.2 


1.5166 


Crll 


0.6±0.3 


5190.3 


2062.7 


1.5166 


Znll 


5883.4 


2344.2 


1.5100 


Fell 


1.4+0.5 


5983.0 


2382.8 


1.5109 


Fell 


2.2±0.5 


6497.5 


2586.7 


1.51 


Fen 


11±2 


6527.0 


2600.2 


1.51 


Fen 


7031.1 


2796.3/2803.5 


1.5112 


Mgn 


5.2+0.3 



